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Calculation of the elastic scattering parameters in an ultra-cold Fermi-Bose and 

Bose-Bose Pr vapor 

M. Kemal Oztiirki-Q and Siileyman OzQcliki-Q 

'Physics Department, University of Gazi, 06500 Teknikokullar, Ankara, Turkiye 

(Dated: Today) 

The calculations of the elastic scattering properties of ultracold-francium isotopes are reported in 
the detail. A parametric model potential for the triplet molecular states of the Fr2 is represented, and 
the scattering lengths at and the effective range r^ are calculated using WKB and Numerov methods 
for Fr-Fr in the triplet state. The convergence of these scattering properties as the depending on 
a Kq parameter and core radius is also investigated using two method as Quantum Defect Theory 
and analytical formula of symskovisky. 

PACS numbers: 32.80.Cy, 34.20.Cf, 34.50.Pi 



I. INTRODUCTION 



In the Bose-Einstein condensation of alkali metal 
atoms, about francium atom has been known veiy little 
knowledge except trapping and cooling until now [j] . Be- 
cause it has too much energy coming from nuclear reac- 
tion and no stable long-lived isotopes. Ultracold atomic 
collisions gained renew an important have came out a ne- 
cessity of research as result of the atomic trapping and 
cooling when energy levels, hyperfine structure intervals 
and isotopes shifts for Fr isotopes is intensely studied in 
detailed^, ^]. Calculations of ultracold elastic scattering 
properties for this atom haven't been seen in literature 
until now. Reason of this is that experimental work on 
Fr has been lacking mostly due to no stable long-lived 
isotope IJ, |3- Collisions of alkali metal atoms except Fr 
atom at a few mill degrees Kelvin temperatures have been 
studied recently in many experiments |a, 0, H, M, U^ Hll ■ 
Theoretical work about scattering on Fr has been scant 
and only a few published papers exist |1^. 

The character of the interaction between atoms in low- 
temperatures is determined by signs (±) and magnitude 
of the scattering length a: for bosonic and fermionic 
atoms with a > Bose condensate is stable, for a < 
it is unstable|13|. The Pauli Exclusion Principle 
(which was determined experimentally) states that no 
two fermion particles can occupy the same state at the 
same time. They must have some way of being distin- 
guished, whether by location, spin state, or some other 
property. That means that if one fermion is in a local 
ground or minimum energy state, the next fermion in 
the area must be in a higher energy state. For bosons, 
however, the Pauli Exclusion Principle is irrelevant by 
definition - so all of the bosons can be in the same state 
at the same time. They don't have to be distinguishable 
from each other. When this happens, a Bose-Einstein 
Condensate is stable. Accurate calculations of s-wave 



scattering length and effective range for diatomic poten- 
tials are important due to the elastic collisions in these 
temperatures dominated by s-wave or low wave scatter- 
ing. The accurate determination of scattering length de- 
pends on choosing the interaction potentials. Collision 
processes at near- zero temperatures (~ fiK) are sensi- 
tive to the details of the interaction potentials between 
the colliding systems. 

In this paper, we have calculated scattering lengths a, 
effective ranges r^ and cross section for collisions of Fr-Fr 
isotopes at ultracold temperatures with the semiclassi- 
cal method (WKB) 14] and also solved the Schrodinger 
equation using Numerov method numerically ll5|]. The 
interaction of francium isotopes in triplet ('^X^m ) ^^ ^^~ 
scribed by short and long range potentials. This inter- 
action potential V{r) for Fr atom was parametered by 
using coefficients found in the literature. These calcu- 
lated scattering properties using this potential are com- 
pared for two methods. We investigated the convergence 
of some scattering properties obtained with help of these 
methods as the dependence of the core radius and a K'^ 
parameter using Quantum Defect Theory|l6j and the an- 
alytic calculations of the scattering lengths obtained by 
Szmytkowski 17]. 



II. PROCESSES 

The scattering length is defined from the asymptotic 
behavior of the solution of the radial Schrodinger equa- 
tion at zero energy: 



dr^ 



2i2V{r) + fc2 



1(1 + 1) 



y/«=o 



(1) 



in atomic unit. Asymptotic behavior of the wave function 
is 
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y(r) = s r + so, 



(2) 



as r ^- oo, where s and sq are constants. The so 
length is given by 



s 

a = , 

So 

here coefficients are obtained by using WKB 

from semiclassical behavior and with help of tl 

solution at zero-energy, eq. 13 is transformed to t 

EF 



m 



a = a{l — tan{iy9 — ^r/S}} 

in atomic unit, where a is the "mean" or "typici 
tering length determined by the asymptotic beh 
the potential through the parameter 7 = y/2jlK 
duced mass, and k = Cg is Van der Waals const 
atom-atom interactions: 



a = v/27r(3/4)/r(l/4) 

Scattering length also depends a semiclassical 
calculated at zero energy from classical turning 
where y{ro) = 0, to infinity, 
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FIG. 1: Changes of scattering length as a function of the 
semiclassical phase. The solid point in the left shows triplet 
scattering length for Fr - Fr isotope and other is the 
result of 1% changes in phase. 



non-zero potential) with the boundary condition UQ(r —f 
00) = as u{r) ~ ua(r) at r —> 00. Then, effective range 
can be written as Il9l 



(p= ^y-2^iV{r)dr. 



It also determines the total number of vibratic 
els with zero orbital angular momentum, Nbd — 
5/8} -|- 1, where { } is the integer part. When 
ference in brackets is just below an integer the 
ing, eq. ^ is anomalously large negative (see fi 
\a\ >> a, which correspond to presence of a virti 
at e = h'^/{2fia^), and when it exceeds an intej 
margin, a is very large positive, due to the existe 
weakly bound state. Unlike 7 and a, the phase : 
depends strongly on the actual shape of the inte 
potential well. When the phase is large, ip/n >; 
scattering length is very sensitive to the slightest 
of the potential 18]. It shows clearly that in figu: 
change of 1% percent in phase gives the less sc 
length as 161.238a.u. instead of 434.477a.u. fo 
2i3pj, jgotope interaction. So, the accurate detern 
of potential is very important for scattering proj 

The effective range r^. can be written in terms 
energy solutions of the partial- wave equation if uo{r) is 
the solution of the partial-wave equation at fc = for 
zero-potential, and normalized as 



»ji. zjci-^y- 



Wo 



sin(fcr -|- i5o) 
sin So 



(7) 



as fc ^ 0, and if uo(r) is the normalized solution of the 
radial equation^for the s partial wave at zero energy (for 



= 2 



"oW -u{r)j 



dr. 



(8) 



This integral converges provided u{r) approaches uo(r) 
rapidly enough as r ~f cx3. This requires V{r) to de- 
crease faster than r^^. For — C6/r'"potential type, this ex- 
pression is adjusted by Flambaum et al. JLSJ using WKB 
method and the exact solution of the radial equation at 
zero-energy to form 



re 



r(i/4) 
r(3/4) 



^ r(3/4) 47 



a r(l/4) a2 



(9) 



Then, the mean effective range may be rewritten by re- 
placing a ^ d into eq. |5| as 



- %/7 



\/2r(l/4) 4^7 r(3/4) 4^/27' 



r(3/4) 



r(i/4) 



a^ 



(10) 



The low-energy scattering is dominated by the contri- 
bution 1 = 0. At values of k close to zero, the I = phase 
shifts (5n can be represented by a power series expansion 
infc [Ifll 



fccot(5o = h -refc^ ' -'''-^ 

a 2 



(fc3) 



(11) 



III. RESULTS AND DISCUSSIONS 



given by 



A. Potential 

The interaction between Fr atoms like the other al- 
kah metal atoms has two features: First, the potential 
at large distances behave as an inverse power of the in- 
teratomic distance, v — —K/r" with n = 6 for spher- 
ically symmetric two atoms in their ground state. The 
asymptotic parameter k = cg is known quite well for most 
atomic pairs of interest. Second, for alkali metal atoms 
other than hydrogen and helium the potential curve is 
generally quite deep and also the electron-exchange part 
of the atomic interaction is repulsive at smaller distance 
as for the singlet and triplet terms of alkali- metal atoms. 
"Deep" here means that the wave functions of the atomic 
pair oscillate many times within the potential well, even 
at very low collision energies. The interatomic potential 
supports a large number of vibration levels. This latter 
feature enables one to use the semiclassical approxima- 
tion to describe the motion of Fr atoms within potential 
well [23. 

Alkali-metal atoms have uncompensated electron spin 
s = 1/2 in the ground state. Therefore, the interaction 
between two alkali atoms in the ground state results in 
the formation of a diatomic molecule, which is described 
by two terms corresponding to the total spin Sab = Sa + Sfe 
of the system equals to 1 or 0. The state with s^fc = 
and 1 corresponds to the singlet state and triplet state, 
respectively. The probability of the spin exchange is de- 
termined by both the dependence of potential curves on 
the internuclear distance r and the splitting between sin- 
glet and triplet terms. The exchange interaction is ana- 
lytically calculated using the surface integral method by 
Smirnov and Chibisov|22J. which yields 



w.M = -K- 



-Vg) = J{A,a,f3)r' 



-lg-(a-(-/3)r^ 



(12) 

Where a = /3 is the ionization energy of francium 
atom, and J are constant. In this equation, J is cal- 
culated from the surface integral given in the Ref . |23| . 
Also we have obtained the value of normalization con- 
stant, A by joining in the range of distances from the 
nucleus the wave function y(r), obtained by the Hartree- 
Fock method, and the asymptotic value of the wave func- 
tion of Eq. (O for an electron far from the nucleus but 
in the Coulomb field of the atomic core. For Fr-Fr, these 
values of constants are given in Table HJ 



^long-range — Zj J c\' J 



(13) 



which is the sum of Van der Waals terms. The C3 and 
Cg Van Der Waals values are given in table |l] had been 
calculated by Derivianko in literature 23]. The potential 
of Fr-|-Fr atom interaction consists of above expressed 
two terms, the sum of Eqs. ll2l and [T!n is given as 



y2.^B+i(^) = V{r)i 



ong — range 



+ (-l)^-T/(r), 



(14) 



Here, for triplet state. Sab is taken as 1 value. This 
potential can be usefull for other studies. 



B. Scattering length and effective range of Fr 

The scattering lengths of francium isotopes calculated 
from Eq21 using WKB method for V{r) are given in the 
Table ^ In order to calculate numerically scattering 
length, we obtained s-wave phase shift Sq using the po- 
tentials in eq^5 then we do it by solving radial equation 
using Numerov algorithm at e = fc^/2^ in atomic unit 
and finding Sq from the asymptotic behavior of the wave 
function y{r) = sin(A;r -|- (5o). The phases at small mo- 
menta k are used to extract the scattering length numer- 
ically from a — \iui(ta,ndo/k). The scattering lengths ob- 

tained via this method are given in the numeric columns 
of Table lUl 



TABLE II: Scattering lengths of the samespecies francium 
isotope interactions in the triplet state. x,y gives isotopes 
cases. 



Triplet scatterin 


g lengths 


at{a.u.) 


IsotopesC'Fr- 


ypr) 


WKB 


Numerov 


•2UH-'2UH 




-52.3151 


-52.303 


209-209 




-134.209 


-134.229 


210-210 




-332.242 


-332.672 


211-211 




-1779.55 


-1779.490 


212-212 




932.576 


932.598 


213-213 




434.477 


434.481 


214-214 




301.464 


301.461 


215-215 




237.152 


151.911 


216-216 




197.235 


197.226 


223-223 




52.0355 


52.035 



TABLE I: Potential coefficients for Francium-Francium atom 
in triplet state. 



A 



C3 Cfi 



0.500 4.480x10"'' 0.5471 4.437 5256 



The long-range part of the interaction potential is 



Some parameters are important in the determination 
of a El for WKB method. We also obtain these param- 
eters, the mean scattering length a and the zero-energy 
semiclassical phases (p, given in Table IIIII using values 
of asymptotic parameter 7. Here, the mean scatter- 
ing length with asymptotic behavior of the potential in 
Ea. (|14|) has been calculated from equation 



TABLE III: Mean scattering length a, semiclassical phase 
shifts (p and numbers of the vibration bound states of same- 
species isotope interactions in the triplet states. 



("Fr- 


-^Fr) 


awKB 


f Nt 


■AUH- 


-2U« 


101.173 


73.6373 22.8145 


209- 


-209 


101.294 


73.8139 22.8707 


210- 


-210 


101.415 


73.9904 22.9269 


211- 


-211 


101.536 


74.1662 22.9828 


212- 


-212 


101.656 


74.3419 23.0387 


213- 


-213 


101.775 


74.5170 23.0945 


214- 


-214 


101.895 


74.6922 23.1503 


215- 


-215 


102.014 


74.8667 23.2058 


216- 


-216 


102.133 


75.0412 23.2613 


223- 


-223 


102.952 


76.2502 23.6462 



As shown in the Table |nl the scattering lengths calcu- 
lated for the collisions of Francium pairs is in excellently 
agreement for two method. 

Even our results shows that the stability of large 
condensates requires repulsive interactions (positive a), 
whereas for attractive interactions (negative a) it is un- 
stable, only a finite number of atoms can be found in 
condensate state in a trap. Due to changes in the V(r), 
scattering length is found as negative value, which lead to 
a condensate triplet state where the number of atoms is 
limited to a small critical value determined by the mag- 
nitude of a . In contrast, we have observed the posi- 
tive scattering lengths that produce stable condensates 
The calculations given in Table ITTl shows that all atoms 
are found in the same energy levels If interaction atom 
species are boson, or they display a adopt behavior due to 
Pauli exclusion principle If interaction atom species are 
fermion. Also, from scattering lengths given as a func- 
tion of the semiclassical phase shift in the eq. Ql , we say 
that the minor changes in phase shifts have strongly af- 
fected the scattering amplitude. It is also remarked that 
for 2i2pj._2i2pj.^ ^t^^ analysis in the direct determination 
of scattering length is weakly dependent on the number 
Nb of bound states supported by ab-initio potential. In 
contrast, mass scalling from one isotope to the others 
depends more strongly on Nf,. 

As shown in table E] scattering lengths are unstable 
against collapse from isotope 208pj,_208pj. ^^ 2iipj,_2iipj, 

Here the both boson-boson and fermion-fermion interac- 
tions are unstable. Same case is valid for the other stable 
isotope interactions. 

Also, the effective ranges calculated using Numerov 
and WKB methods in the function of phase for the V{r) 
potentials are presented in Table IIVI The close agree- 
ment between the calculations of r^ from Eqs.JSJ and 
(|ll|l confirms the accuracy of the numerical integration 
of the partial-wave equation. The size of the scattering 
lengths and effective ranges is closely related to the po- 
sition of the last vibration bound states of the energy 
curves, as can be anticipation by inspection of Eq.© 
and number of vibrational levels ns(<p/7r) which, consis- 
tent with Levinson's theorem, show that as the binding 



energy of the highest level tends to zero, the scattering 
length tends to ± infinity|2J|. 



TABLE IV: Effective ranges of the samespecies francium iso- 
tope interactions in the triplet state. x,y gives isotopes cases. 



Effective ranges rt 


(^Fr.apr) 


rwKB rNumerov 


WH- 


-208 


3645.27 3635.894 


209- 


-209 


1078.46 1078.008 


210- 


-210 


531.721 531.615 


211- 


-211 


332.01 332.150 


212- 


-212 


239.005 239.243 


213- 


-213 


190.433 190.127 


214- 


-214 


164.266 164.168 


215- 


-215 


151.738 151.738 


216- 


-216 


149.197 149.197 


223- 


-223 


1463.56 1463.560 



C. Convergenge of scattering length and effective 
range 

Analytical calculations of scattering lengths are im- 
portant in investigation of convergence of the scattering 
length and the effective range. In many approaches used 
to solve collision problems in atomic physics, the tree di- 
mensional configuration space is divided into two regions 
separated by a spherical shell (a core boundary) of radius 
p. In the inner region (r < p) the short-range interaction 
between two colliding particles is very complicated and 
a scattering equation must be solved independently for 
each combination of particles. In contrast, it p is chosen 
sufficiently large the scattering problems in the outer re- 
gion (r > p) may be reduced to potential scattering with 
the long-range potential accurately approximated by sim- 
ple analytical expression. A numerical solution in this 
region is usually easily approachable 2 5j . Two of these 
methods are analytical calculations derived by Radoslaw 
Szmytkowski for inverse problems '25l and Quantum De- 
fect Theory described by Gao 16]. 

Radoslaw Szmytkowski solved the radial schrodinger 
equation at zero energy for inverse type of long-range in- 
teraction analytically with solutions expressed in terms 
of the Bessel, Whittaker and Legendre functions, respec- 
tively. They derive exact analytical converge formulae 
for the scattering lengths. The expressions depend on the 
short-range scattering length, core radius and parameters 
characterizing the long-range part of the interaction. In 
Table We obtained the convergence of the ultra cold 
scattering length of Francium atom using the formulae of 
Szmytkowski. These formulae are given in that paper as 
eqs. (75), (73). As shown in the table, this method con- 
firms convergence of the scattering length for ^^^Fr-^^^Fr. 
Here, only, the convergence of the scattering length calcu- 
lated from WKB method is given, and the same behavior 
can also be seen in the numerical calculation. 



TABLE V; Convergence of scattering length calculated by 
WKB with the formulae of Szmytkowski for ^^^Fr-^^^Fr. 



p 


ao 


102 


10 


-73.2943 


-6.19E+09 


100 


53.1953 


-46796 


500 


929.775 


929.775 


1000 


932.426 


932.426 


5000 


932.569 


932.569 


10000 


932.575 


932.575 


50000 


932.576 


932.576 


Infinity 


932.576 


932.576 
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Quantum defect theory of atomic collisions is presented 
by Gao. Based on the exact solutions of the Schrodinger 
equation for an attractive l/r^ potential, the theory 
provides a systematic interpretation of molecular bound 
states and atom-atom scattering properties and estab- 
lishes the relationship between them. He finds a defini- 
tion for the scattering length and the effective range and 
s wave at zero energy as 



ai=o 



ri=o 



[r(l74)]^ -pn-P6 



-P6 



7(yp— Pe 



3ir 



\K» 



(15) 



where (3q = (2iJ,Ce) and X'^is the analytic function of 
energy. 

Figure [21 shows the scattering length and the effective 
range curve plotted using the eqs. ((T3|l and ones (dotted 
points) obtained by the mean WKB formulas. 

Our graphical representations show that as K^ goes 
infinity, the QDT results converges excellently to values 
obtained by WKB method in Fig. [21 Even, it may be 
said that the different calculation of these properties gives 
similar results for any ultracold atom-atom collision in 
these two method. 
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FIG. 2: Convergence of the scattering length and the effective 
range calculated with QDT as a function of K'^ parameter 
when K'^ goes infinity. 



multiplied with 2 value due to symmetric wave function 
If systems are bosonic one. Table fVll shows zero energy- 
cross sections are adopted for two method. It concludes 
that zero energy cross sections is expected to be smaller 
than that of other alkali metal atom, because of it is a 
heavier atom 



D. The zero energy cross sections 

It is also known in more tex book that the low energy 
scattering is always dominated by the I = partial wave 
the corresponding phase shift being given by 



Sq = ka (16) 

The triplet elastic cross sections can be defined bv[2^ 



47r ^ 

1=0 



2 xt 



(21 + 1) sin^ S, 



(17) 



where / — 0, sin^o = Sq in limit of zero-energy and Eq. 
nn is transformed using eq. El to cr*; = Atto? . This 
result is valid for fermion-fermion interaction, but it is 



TABLE VI: Triplet cross section of different francium isotopes 
calculated for two method. The numbers in brackets denote 
multiplicative powers. 



Triplet cross section(»n 



(^Fr-^Fr) WKB 



Numerov 



-'Am 

-209 



209- 

210-210 

211-211 



215-215 
216-216 



1.92617(-16) 1.92528(-16) 

6.33833(-16) 6.34022(-16) 

7.76874(-15) 7.78886(-15) 

1.11438(-13) 1.11430(-13) 

6.12083(-14) 6.12111(-14) 

6.64271(-15) 6.64283(-15) 

6.39606(-15) 6.39593(-15) 

1.97909(-15) 8.12064(-16) 

2.73785(-15) 2.73760(-15) 

9.52820(-17) 9.52802(-17) 



IV. CONCLUSIONS 

The elastic scattering properties for the collision of the 
francimn atoms in a limited range of moment (fc < a~^ ss 
O.Ola.M. for francium-francium isotopes) at low tempera- 
tures are sensitive to the details of interaction potentials. 
Scattering properties such as scattering length, effective 
range and cross section have been computed using semi- 
classical and a numeric methods for the V{r) potentials 
as dependence of cutoff radius adjusted by comparing 



with the experimental potential for ultra-cold francium- 
francium isotopes collision. We investigated the conver- 
gence of these scattering properties as the depending of 
core radius and K'^ parameter using Quantum Defect 
Theory and analytical calculations derived by Radoslaw 
Szmytkowski for inverse problems. Cross section was ob- 
tained as ~ 10^ ^^m^ at low energy. The phase shifts, an 
intermission parameter for scattering length and effective 
range, has linear manner at small momenta k < a~^. 
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